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1. Introduction 

The Gross-Pitaevskii equation (GPE) is one of the basic models in the theory of Bose- 
Einstein condensate (BEC)[1]. In dimensionless form this equation reads 

i% = A^-\/(x)^ + (j^|^| 2 , (1) 

In BEC context GPE describes an atomic cloud in the mean-field approximation. The 
function \I/(t,x) corresponds to macroscopic wave function, A is the Laplace operator 
and V(x) is a trap potential for the condensate to be confined. The value of a governs the 
sign of nonlinearity: the case a = 1 corresponds to attractive interatomic interactions 



and a — — 1 corresponds to repulsive interatomic interactions. Both the cases a — ±1 
are of interest for physical applications [2, 3]. Eq.(l) conserves the quantities 

E=^J (|W| 2 - V(x)|^| 2 -||^| 4 ) rfx (2) 

N = J \m\ 2 dx (3) 

which have the sense of Hamiltonian and number of particles correspondingly. An 
important class of solutions of GPE are stationary nonlinear modes defined as 

tt(f,x) = e^V(x) (4) 

ip(x) ^ as |x| — > oo (5) 

In this case satisfies the equation 

Aip + (tu- V(x))^ + cn/#| 2 = 0, (6) 



The parameter u in terms of BEC corresponds to the chemical potential. The ground 
state solutions of (5)-(6) (i.e. node-free solutions which minimize energy functional 
for Eq.(l)) are of primary importance for BEC applications [4, 6, 5]. Apart from 
them some non-ground (i.e. excited) nonlinear modes also have been reported (see 
e.g. [7, 9, 11, 12, 13]) as well as physical mechanisms for their generation [8]. Except 
some specific examples of the potential V(x) when exact solutions of (5)- (6) are available 
(see e.g. [14]) in all the cases the nonlinear states have been found numerically. 

Typically, numerical calculation of nonlinear modes needs some initial guess for 
subsequent iterative procedure. In order to get a priori information about the shape 
of nonlinear states various approaches have been used. In [9] the branches of nonlinear 
modes for the case V(x) = |x| 2 have been found by continuation from the limit i/j — > 
where the nonlinear term can be neglected 

At/j + (u - |x| 2 )V = 0. (7) 

Therefore, the shapes of these states can be described by the set of Gauss-Hermite 
functions. In ID case and for a — 1 it has been proved rigorously [10] that there are 
infinitely many branches of nonlinear modes which bifurcate from zero solution of (6) 
at discrete infinite set of values of uo corresponding to eigenvalues of linear problem 
(7). Following [11, 12], it is quite natural to call these modes the modes with linear 
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counterpart. It has been shown in [11, 12] that for some potentials V(x) nonlinear modes 
without linear counterpart also can exist. In the case of multi-well potentials these states 
have been found by continuation from asymptotical limit of large nonlinearity [11, 12]. 

However, the following question remains: can one list all nonlinear states 
corresponding to given potential V(x) with oj lying within some range? To best of 
our knowledge, up to the moment there is no answer to this question even in the 
most explored ID-case of quadratic potential V(x) = x 2 . The aim of this paper is 
to make a step toward answering this question. Specifically, for a = —1 we explain 
how to arrange numerical calculation in such a way that all nonlinear modes described 
by Eq.(6) for u lying within a given range would be found. In the case a = 1 the 
method gives a picture which is complete, as we believe, but this fact has not been 
proved. We illustrate our approach on ID case of Eq. (6) with quadratic and double- well 
potentials (the applications to radial 2D and 3D cases are straightforward). The base 
for this approach is given by exact statement on one-to-one correspondence between the 
solutions of Eq.(6) which tends to zero as x — > +oo and linear coefficient C + which 
determines their asymptotics in this limit. 

So, in this paper we consider ID version of Eq.(6), specifically 

ijj" + V(x))iP + aip 3 = (8) 

where prime denotes differentiation with respect to x. We assume that the potential 
V(x) is even, V(x) = V(— x) and grows when \x\ — > oo. The attention is focused on two 
particular examples of the potential, 

V(x) = x 2 , (quadratic potential) (9) 
V(x) = ax 4 + (3x 2 , a > 0, (3 < (double-well potential) (10) 

The solutions of Eq.(8) are supposed to be real and satisfy the localization condition 

lim %l)(x) = (11) 

rr^ztoo 

The paper is organized as follows. In Section 2 we describe the method which 
we use. In Section 3 we apply this method to Eq.(8) with quadratic and double- well 
potentials. Section 4 contains summary and discussion. Some auxiliary statements and 
technical details can be found in Appendices A-C. 



2. The description of the approach 

2.1. The idea of the approach 

The approach which we propose is an extension of one used in [6, 20] for numerical 
construction of nonlinear modes for Eq.(6). Let us denote by S + the set of all solutions 
of Eq.(8) which satisfy the condition 

lim ip(x) = (12) 
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When x — > +00 the nonlinear term in Eq.(8) can be neglected and the asymptotics of 
solution from S + can be found from linearized equation 

ijj" + V(x))ip = (13) 

The theory of Eq.(13) (see e.g. [15, 16]) says that for a wide class of potentials V(x) it 
has unique (up to constant factor) solution i]){x) such that i]){x) — > when x — > +00. If 
V(x) > uj for x > x for some Xo large enough (and fixed) then the asymptotics of this 
specific solution is described by the formula [16] 

-J y/V{t)-W dt 

^) = \'l 7rT ^ (l + o(l)). *-+oo (14) 

y/V(x) - UJ 

Then the solution of Eq.(8), ip(x), which satisfies (12) has the asymptotics 



- / y/V(t)-W dt 

^{x) = € 'L^^ (C+ + o(l)), x^+00 (15) 

Here C + is some constant. The next step is to assume that there is one-to-one 
correspondence between the solutions from S + and the constants C + G R which comes 
from formula (15) (this fact will be rigorously justified below). This implies that the 
solutions S + constitute one parameter set with free parameter C + G R. Let us denote 
a solution of Eq.(8) which satisfies the asymptotical relation (15) by ip + (x,C + ). By 
construction, the solutions ip+(x; C + ) tend to zero as x — > +00; however, generically, 
they do not tend to zero as x — > —00. 

The same reasoning holds for the solutions of Eq.(8) which tend to zero as x — > —00 
(the set of these solutions will be denoted by S'_).The asymptotics of these solutions is 
of the form 

f°y/v(t)-u> dt 

^ x ) = ^mn= ' ( C7 - + ( 1 ))' ( 16 ) 

i/V[x) - uj 

where C_ is some constant. Assume that S- is one-parameter set with free parameter 
C-. Let us denote a solution of Eq.(8) which satisfies the asymptotical relation (16) by 
ip-(x, C-). Evidently, a localized solution ip(x) which satisfies (11) must belong to both 
sets, S + and S-, and in the point x = the relations must be valid 

^(0) = V-(0,C_) =^ + (0,C + ); V'(0) = ^(0,C_) =^;(0,C+); (17) 

In fact, (17) is a system of equations to determine unknown constants C_ and C + and, 
consequently, a solution ip(x). In practice, the values ip + (0;C + ) and ■?/ ; + (0>C+) can be 
calculated numerically (using, for example, Runge-Kutta method) by solving the initial 
value problem for Eq. (8) starting from some large enough value x = x^ where the values 
■ifj+i^Xoo] C + ), ^+(^00; C+) can be approximated by means of asymptotical formula (15). 
In the same way the values t/>_(0; C_), ^-(O; C_) can be calculated. 
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2.2. Mathematical basis 

The heuristical arguments given in paragraph 2.1 can be justified by exact statement 
which describes the set of solutions of Eq.(8) which tend to zero as x — > +oo. This 
can be done using "asymptotic integration" of Eq.(8) (see e.g. [19]). The following 
statement is valid: 

Theorem L Let there exists some x such that 

(i) V(x) — uj > e > for x > x ; 

(ii) V{x) G C 2 (x ,+oo); 
(Hi) f°° n )Y, ^ttto dx < oo 

» 7 Jxo (V(x) — w)- i >' ! 

Let ipi (x) be a fixed solution of linear equation (13) such that ipi (x) — > as x — > +oo 
which has asymptotics (14). Then 

(a) for any solution ip(x) G S + there exists a limit 

lim ^^ = C+ <oo (18) 

(b) for any C+ G R there exists an unique solution ip+(x; C+) G S + of Eq.(8) such 

that 

|im M^c ± ) =c+ 

ipi{x) 



(c) if ip+(x; C + ) exists in some neighbourhood U of a point (x; C + ) then in U there 
also exists partial derivative ^+^ c +) which is continuous with respect to C + . 

Theorem 1 establishes one-to-one correspondence between S + and R; specifically, 
it says that the set S + can be parametrized by real parameter C + G R. The point (a) 
of Theorem 1 has been used for the search of localized solutions by means of shooting 
method (see e.g. [20]). We make use of the point (b) to construct general picture of 
branches of solutions of Eq.(8). The point (c) will be used for proving Lemma 1 in 
Appendix C. 

Theorem 1 reflects the well-known relation between the behavior of linear dynamical 
system and its nonlinear perturbation in a vicinity of its zero solution. In spite of bulk 
of literature devoted to this subject (see, for instance, [17] and references therein) we 
failed to find this statement exactly in the form we need. Therefore we adduce a sketch 
of proof of Theorem 1 in Appendix B in a little more general form. Evidently, analogous 
statement can be made for the solutions of Eq.(8) which tend to zero as x — > — oo. 

2. 3. Practical implementation of the method and visualization of the results 

2.3.1. (ip(0),tp'(0)) -pictures. Let the potential V(x) be even, V(x) = V(—x). Let 
uj be fixed and x = Xoo is some value large enough where the solution of (8) can be 
approximated well by the asymptotics (15). Assume that the numerical continuation 
of the solutions ip+(x; C+) from x = x^ to x = can be fulfilled for each C + in some 
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interval C + G [— C + ; C + \ (using, for instance, Runge-Kutta method). The result can be 
represented as a segment of curve 

7+ : {(^ + (0;C + );^(0;C + )),C + e [-C + ;C + ]} 

on the plane Since C + = corresponds to zero solution, this curve passes 

through the point (0; 0) of this plane. Moreover, due to the fact that the equation 
(8) admits the symmetry ip — > —ip, this segment of the curve is centrosymmetric with 
respect to the origin. In practice this means that the calculation should be fulfilled for 
[0; C + ] only. In analogous way, let us consider the segment of curve 

7- : {(V>-(0;C_);^(0;C_)),C_ G [-C.;C.]} 

constructed by means of numerical continuation of the solution ■0_(x;C_) from the 
value x = — until the value x — 0. Since equation (8) is invariant with respect to the 
change x — > —x, if C_ = C + the curve 7_ is symmetric about the ip' axis to the curve 
7 + . Then, each point of intersection of 7_ and 7 + corresponds to a solution of system 
(17) and, consequently, to a solution of Eq.(8) with the condition (11) (see Fig.l). 

Due to the symmetries of the curves 7_ and 7 + , the points of intersection of 7+ 
with the axes ip and tp' are also points of intersection of 7+ with the curve 7_. These 
points correspond to specific localized solutions of Eq.(8); specifically, the points of 
intersection of 7+ with the axis ip correspond to even localized solutions of Eq.(8) 
whereas the points of intersection of 7+ with the axis ip' correspond to odd localized 
solutions of this equation. The points of intersection of 7_ and 7+ which do not belong 
to the axes correspond to non-symmetric localized solutions of Eq.(8). 

2.3.2. (C+, uo)- diagrams. Let us consider now branches of localized solutions of Eq.(8) 
which appear when the parameter u varies. Due to Theorem 1, if 00 is fixed, each of the 
localized solutions of Eq.(8) can be unambiguously defined by the value C + determined 
by its asymptotics. Therefore, it is convenient to represent the branches of localized 
solutions on the plane (C+,u) marking for each uj those values C + for which localized 
solutions exist. In such a way any branch of solutions corresponds to a curve Y on 
the plane (C+,u). An advantage of this graphical representation (in comparison, for 
instance, with more traditional (N,u) representation where N is given by (3)) is that 
each point of the plane (C + ,ou) corresponds to at most one localized solution of Eq.(8). 
This means, for instance, that if two curves, corresponding to two branches have a 
common point, then there exists a localized solution which belongs to both branches. 
This is not the case for (N, a>)- diagrams. 

Consider now any segment of curve T which has no intersection with the axis 
C+ = 0. Simple reasoning shows that the number of zeros of solutions ip(x) is constant 
along this segment of T. Indeed, consider a continuous deformation of solution tp(x) 
which occurs when passing along this segment. A zero x = x of ip(x) can disappear 
by (i) merging with one or several neighboring zeros; (ii) by merging with x = ±00 
("goes to infinity"). The case (i) implies that the branch includes a solution with 
ip(xo) = ip'(xo) = which is zero solution. Consider the case (ii) i.e. assume that zero 
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x = Xq merges with x = +00. Then for some solution of the branch the parameter C + 
is zero which implies that the solution is identically zero. Similar arguments allow to 
discard also the case when (i) and (ii) hold simultaneously. 

2.4- Some additional statements for the case of repulsive nonlinearity 

In the case of repulsive nonlinearity Eq. (8) reads 

tP" + V(x))ip -ip 3 = (20) 

It turns out that "most part" of solutions of this equation are singular i.e. they tend to 
±00 as x approaches some finite value. This fact allows to fulfill more complete analysis 
of localized solutions for Eq.(20). It can be based on the following statement: 

Lemma 1 Let there exists a solution ipo(x) = i/)+(x, C°) G S + , C° > of equation 
(20) such that: 

(i) ipo(x) is defined on [xq\ +00) where x > and 

lim ipo(x) = +00 

(ii) ipo(x) > for x G (x ; +00); 

(Hi) for x G (x ; +00) the following inequality holds: 

uj-V(x) -3ipt(x) <0 

Then any solution ip + (x,C + ) G S + for C + > C° has a singularity at some point 
x = x(C+) > x , \im x ^ x{c+)+0 ip + (x,C + ) = +00. 

The proof of Lemma 1 is postponed in Appendix C. Lemma 1 allows to stop the 
scanning C + when C° satisfying the conditions of Lemma 1 is achieved. The window 
of scanning in this case is finite, so one can conclude that all the solutions which exist 
for given uu have been found. However, for practical implementation of Lemma 1 in 
numerics one needs a condition to assure that the solution which has been found is 
singular. 

Let us introduce the value 

fl = sup(c<j — V(x)) 

and assume that |Q| < 00. If SI < then one can prove that Eq.(20) has no localized 
solutions. If SI > the following statement holds 

Lemma 2. Let Q > 0. If a solution ip(x) of Eq.(20) satisfies at some point x = Xo 
the condition \i/j(x )\ > VTi then this solution is singular. 

The proof of Lemma 2 is also postponed in Appendix C. Lemma 2 allows to stop 
the shooting procedure once the value ip = y/Ti is achieved. 
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Figure 1. (ip(0), V>'(0))-pictures for Eq.(8) with V(x) = x 2 . Left panel: repulsive 
nonlincarity, ui = 4; right panel: attractive nonlinearity, cu = 0. 



3. Numerical results. 

S.i. The quadratic potential. 

3.1.1. Repulsive nonlinearity First, we applied the approach described above for the 
equation 

tp" + (uj - x 2 )i> - i> z = (21) 

taking < uo < 12. According to the formula (14) the asymptotics of the solution ip(x) 
as x — > +oo is 

^(x) ~ C+x^^-^e-T (22) 

The starting point for the shooting, taken equal to 10. The results were re- 

calculated for Xqo = 12 and no significant change was observed. The zoom of scanning 
with respect to C± was from C± = until the value C±^ where the conditions of Lemma 
1 held (the value C±^ has been found for all considered uo). This guarantees that all 
localized solutions within this range of uj were found. 

The Fig.l, left panel, shows the (^(0), ip'(0))- picture for cu — 4. In Fig.2(a)-(b) 
six branches of nonlinear modes for Eq.(8) with repulsive nonlinearity are depicted. 
All these branches, T™, are originated by localized modes of linear harmonic oscillator, 
corresponding to eigenvalues u = 2n — 1, n = 1 -r 6. The linear limit can be achieved 
by passing C + — > (or, correspondingly, N — > 0). The branches of even and odd 
solutions alternate. It follows from these figures that there are no other branches of 
localized modes within this range of u. The analysis of linear stability for the branches 
shows that the solutions corresponding to the branches Tl and T 2 r are stable whereas 
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Figure 2. The branches of nonlinear modes for Eq.(8) with V(x) = x 2 . 
(a): (C + , w)-diagrams, repulsive nonlinearity; (b): (N, w)-diagrams, repulsive 
nonlinearity; (c) (C+, w)-diagrams, attractive nonlinearity; (d) (N+, c<j)-diagrams, 
attractive nonlinearity. All the modes are originated by the modes of linear harmonic 
oscillator. 



the solutions of the other two branches (and, probably, of all the higher branches) are 
unstable. This completely agrees with results of [21], (see Fig. 1 of [21]). 

3.1.2. Attractive nonlinearity For the case of attractive nonlinearity, 

< + (cj- x 2 )^ + V 3 = (23) 

it is proved [10] that each localized mode of harmonic oscillator, corresponding to 
u; — 2n — 1, n = 1, 2,... originates a branch of nonlinear modes, T™. The asymptotics 
of the solution ip(x) asi^ +oo is described by the formula (22). We study Eq.(23) for 
— 5 < cu < 5. Since in this case Lemma 1 cannot be applied the scanning with respect 
to C± was fulfilled within the fixed range < C± < 10. Fig.l (right panel) shows 
(■0(0), tp'(0))- picture for u = 0. Fig.2(c) represents (C+, w)-diagram and Fig. 2(d) shows 
the corresponding (N, w)-diagram. The four branches r™, n = 1,2,3,4 originated by 
the localized modes of harmonic oscillator have been found and no other modes have 
appeared within the considered range of C± and uj. We suppose that, in general, the 
branches of solutions generated by the modes of linear oscillator constitute all nonlinear 
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localized modes that can be described by Eq.(23). 

Linear stability analysis shows that the solutions of the two lowest branches T\ 
and are stable and the solutions of the branch (and, probably of all the higher 
branches r™, n > 4) are unstable. This conforms with the results of [21], (see Fig.l of 
[21]). At the same time we observed that the two-node solutions of the branch r 3 are 
unstable in small-amplitude limit (for uj* < uj < 5, uj* m 3.83) but they become stable 
for larger amplitudes (uj < uj*). Probably this is the reason of troubles with generation 
of this states by adiabatic variation of scattering length reported in [21]. 

3.2. The double- well potential 

3.2.1. Repulsive nonlinearity . In this case the Eq.(8) reads 

tp" + (u- ax 4 - (3x 2 )i> - i; 3 = 0, a > 0, (3 < (24) 
According to the formula (14) the asymptotics of the solution ip(x) as x — > +oo is 

„ 9± . ex p {-y^x 3 - JL X -J-(u 2 + x- 1 } (25) 
^ y ' x I 3 2^ 4a J J K J 

All the calculations were done for the parameters of potential a — 10, f3 — —20 for 
—5 < uj < 10. The range of scanning for C± was from C± = until the value C±^ 
where the conditions of Lemma 1 held. The value was found for all uj which we 
considered. This implies that all nonlinear modes within this range of uj were found. 

Fig. 3 represents (^(0), ip'(0))- picture for uj = —3. Fig.4(a)-(b) shows four branches 
of localized modes T™, n — 1, 2, 3, 4 which depart from the eigenvalues of linear problem 
(in this case no analytical form of linear modes is available). In addition, the branch 
undergoes secondary pitchfork bifurcation and two branches Tf of nonsymmetric 
solutions arise. These solutions are localized in one of the wells of the potential V(x). 
The linear stability analysis shows that the modes corresponding to the branches T^, 
T^, as well as Tf are stable. The nonlinear modes corresponding to the branch Y 2 . 
are stable until the point of pitchfork bifurcation. After the point of bifurcation they 
become unstable. This is quite expectable due to the results of the papers [25, 26, 27] 
and refines the statement of [12] where it has been suggested that all nonlinear modes 
in the case of repulsive nonlinearity are stable. 

3.2.2. Attractive nonlinearity. In this case Eq.(8) reads 

i>" + {uj- ax 4 - (3x 2 )i) + V 3 = 0, a > 0, (3 < (26) 

The most part of calculations were done for a — 10, f3 — —20. Fig.4 represents all 
the branches of localized modes which were found within the range < C± < 0.4 
and —20 < uj < —4. There are two branches which correspond to modes with linear 
counterpart, and T 2 originated from 1-st and 2-nd eigenvalues of linear problem. At 
some value of uj the branch undergoes pitchfork bifurcation and two more branches of 
non-symmetric solution emerge. Apart from them there is a branch T° of nonlinear 
modes localized on the maximum of the potential situated between the two wells. These 
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Figure 3. (^(0), ■0'(O))-pictures for Eq.(8) with V(x) = 10x 4 — 20x 2 and repulsive 
nonlincarity, ui = — 3. 

solutions were reported also in [12]. When widening the window of scanning with respect 
to C± and u much more branches appear; as a result, the complete picture becomes 
quite complex. The linear stability analysis shows that the solutions of the branch Y 2 a as 
well as the branches are stable, whereas the solutions of the branch T° are unstable. 
As it was in repulsive case, the nonlinear modes corresponding to the branch r* are 
stable until the point of pitchfork bifurcation and are unstable after this point. The 
part of solutions of the branch shown in Fig.4(c)-(d) are unstable. 

The presence of pitchfork bifurcation reflects the fact that Eq.(26) with even 
potential is invariant with respect to transformation x — > —x. The two nonsymmetric 
states ^^(x) G T+ and ip(~\x) G T~ which bifurcate from the principal state 
ipi(x) G are connected by the symmetry -?// + )(— x) = ip(~\x). Adding small 
perturbation which breaks evenness of the potential results in replacing the pitchfork 
bifurcation by saddle- node one (see Fig. 5). Similar phenomenon was reported in the 
paper [26] for the case of the potential composed by harmonic and periodic parts. 

One more interesting phenomenon takes place when changing the parameters of 
the potential. It was found that if a is fixed there exists a critical value of (3 when the 
branches T\ and T° a touch. When f3 passes through this critical value a reconnection of 
branches T\ and T° a occurs (see Fig. 6). This means that by small variation of parameter 
(3 (keeping a fixed) in vicinity of this critical point, a solution with linear counterpart 
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Figure 4. Panels (a)-(b): The branches of nonlinear modes for Eq.(8) with a = 10, 
(3 = —20 and repulsive nonlinearity. (a): (C + , w)-diagram, (b): (N, u)-diagram. The 
branches T™, n = 1,2,3,4 arc originated by the modes of linear eigenvalue problem, 
the branches T^r appear as a result of secondary branching of (the bifurcation is of 
pitchfork type). In (N, u)-diagram the branches T+ and T~ coincide. 
Panels (c)-(d): The branches of nonlinear modes for Eq.(8) with a = 10, (3 = —20 
and attractive nonlinearity. (c): (C+, w)-diagram, (d): (N, u)-diagram. The branches 

and are originated by the modes of linear eigenvalue problem, the branches 
appear as a result of secondary branching of T\ (the bifurcation is of pitchfork type). 
In (N, o>)-diagram the branches T+ and T~ coincide. 

can be transformed into a solution without linear counterpart and vice versa. 
4. Conclusion 

In the paper a method to study nonlinear modes for the Gross-Pitaevskii equation is 
described. It is based on exact statement on one-to-one correspondence between the 
solutions of Eq.(8) which vanish at x — > +oo and a linear coefficient C + which describe 
its asymptotic in this limit. This method can be regarded as an extension of shooting 
method used for calculation of single nonlinear mode (see e.g. [6, 20]) aimed to get a 
complete picture of coexisting nonlinear modes. Combined with some additional exact 
statements which have been proved in the case of repulsive nonlinearity the method 
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Figure 5. (C+, cj)-diagrams for Eq.(8), attractive case, V(x) — 10x 4 — 20x 2 (left 
panel) and V(x) — 10x 4 — 20x 2 + O.Olx 3 (right panel). Breaking of the symmetry of 
the potential results in replacing of pitchfork bifurcation by saddle-node one. 




Figure 6. Transformation of (C+, w)-diagram when the parameters of the potential 
vary. Left panel: a = 10, /? = —20, central panel: a = 10, /? = —15.91, right panel: 
a = 10, /3 = —15.5. Only branches of even solutions are depicted. 

allows to fulfill demonstrative computation to guarantee that all the nonlinear modes 
corresponding to given value of u (chemical potential in terms of BEC) are found. 

Conforming to the approach which we describe we offer a representation of the 
branches of solutions on the plane (C + ,uj) instead of more traditional representation 
on the plane (N,uj), N = J^^ix) dx. The advantage of this representation is that 
each point of the plane (C+, uj) corresponds to at most one nonlinear mode whereas 
the point on the plane (N, uj) can correspond to many nonlinear modes which share the 
same number of particles. Therefore the intersection (or touching) of the curves which 
correspond to the branches of nonlinear modes in the plane (C + ,u) means that some 
bifurcation occurs. This is not the case for (N, uj) representation. The (C + , a;)-diagrams 
calculated in the case of harmonic potential V(x) = x 2 indicate that all the nonlinear 
modes have linear counterpart. In the case of double- well potential the situation is richer, 
due to modes without linear counterpart and secondary branching of non-symmetric 
modes. 

In the paper the approach is applied to ID GPE with potential having one or two 
minima. However it seems to be an interesting problem to apply this method for the 
study of ID GPE with potentials with three or more minima (see e.g. [28]). Also it may 
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be used for GPE with periodic potential (see e.g. [29, 30]) as well as for the potential 
generated by magnetic fields and optical lattice [26]. Applications to quintic or cubic- 
quintic GPE model are also possible [31, 32]. With minor modifications this approach 
can be applied to 2D and 3D radial versions of GPE. 
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Appendix A. Summary of results for the equation ip" — Q(x)ip = 0. 

The equation 

ijj" - Q(x)ip = (A.l) 

(the prime denotes the derivative with respect to x) has been discussed in numerous 
books and papers [15, 16, 18]. Let us set forth below some results about Eq. (A.l) 
which have been used in this paper. 

Statement A.l: (see [15], chapter 2) Let there exists some x and e such that 

Q(x) >e>0 for x > x Q . (A. 2) 

Then a pair of linearly independent solutions of (A.l), ip\(x) and ip 2 (x), can be 
chosen in such a way that for some value a > x the following assertions hold: 

(Al) the functions ipi(x), ip[(x) — > when x — > +oo; 

(A2) the functions ipi(x), ip[(x) have no zeros and are monotonic for x > a; 
(A3) for x > a 

\Mx)\<A 1 e~^ x (A.3) 

for some A\. 

(A4) the solution ip\(x) is unique (up to constant factor) . 
(Bl) the functions ip 2 (x) , ip' 2 (x) — > oo when x — > +oo; 

(B2) the functions ip2 (x) , ip' 2 ( x ) have no zeros and are monotonic for x > a; 
(B3) forx>a 

\Mx)\>Me^ £X (A.4) 

for some A 2 . 

Comment: The points (Al)-(A3) and the fact that the sign of if)"(x) > coincides 
with the sign of ip(x) for [a; +oo) imply that 

ip^x), ^[(x), ip'((x) e L p [a; +oo); p>0 (A.5) 
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Imposing more restrictions to Q(x) results in more precise description for the 
solutions ip\(x) and ^{x) growth or decay: 

Statement A.2: (see [22], Chapter IV, Theorem 14) Let the function Q(x) satisfies 
the condition (A. 2). Let, in addition, the following conditions hold: 

(a) Q(x) G C 2 (a;o,+oo); 

o) 1: dx < 00, 

Then the solutions ip\(x), ^(x) (which exist due to Statement A.l) obey the 
following asymptotic formulas: 

„—S(xo,x) 

Mx)= gv^y (1 + o(1)) ' (A - 6) 

X 

e S(x ,x) r 

^ (x) = gv^ (1 + o(1)) ' s ^ x ) = jVW) dt ( A - 7 ) 

x 

Moreover, the asymptotic formulas (A.7)-(A.6) can be differentiated at least once. 
Comment: Due to (A. 2) S(x , x) — > +00 as x — > +00. 

Comment: Evidently, there is a freedom in choosing the lower bound of integration 
in (A.7)-(A.6). Its change from x = xq to x = x\ results in rescaling of ipi(x) an d i>iix) 

i) X {x) — > e s(xo ' Xl Vi(a:); M*) — ► e~ s{x °' Xl) i; 2 (x) 
Appendix B. Proof of Theorem from Section 2.2. 

The aim of this Appendix is to prove the following statement for the equation 

ip" - Q(x)i> + F(^) = (B.l) 
Theorem L Let there exists some x such that 
(Ql) Q(x) > e > for x > x ; 
(Q2) Q(x) G C 2 {x ,+oo); 
(Q3)!Z^dx<oo 
and there exist some 5 such that 
(NL1) F(t) G C 3 [-5;5]. 
(NL2) F(0) = F'(0) = 0; 

Let ipi(x) be the solution of linear equation ip" — Q(x)ip = such that tpi(x) — > 
as x — > +00 (see Appendix A). Let S + be a set of all solutions of Eq.(B.l) such that 
i/j(x) — > when x — > +00. Then 

(a) for any solution if)(x) G S + there exists a limit 

lim ^\ = C<oo (B.2) 
x^+00 ip-]_[x) 
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(b) for any CgR there exists an unique solution i/)+(x; C) G S + of Eq.(B.l) such 

that 

lim = C (B.3) 

x^+oo ipi{X) 

(c) if ip + (x; C + ) exists in some neighbourhood U of a point (x; C + ) then in U there 
also exists partial derivative dV'+O^cv) w hich is continuous with respect to C+. 

Comment : If for a function F(t) the condition (NL2) holds and F(t) G C 2 [— S;S] 
for some S > then f(t) = F(t)/t has no singularities in \t\ < S and there exists A > 
and Si such that 

\f(t)\ < A \t\, \t\<Si, (B.4) 

The proof of Theorem 1 can be sketched in a form of a sequence of Lemmas. 
First Lemma states that small perturbations of potential in Eq.(A.l) do not change the 
asymptotic (A. 6) of the solution which tends to zero as x — > +oo. 

Lemma B.l Let Q(x) satisfy the conditions (Q1)-(Q3) of Theorem 1 for x > x 
and ipi(x) be a solution of equation tp" — Q(x)ip = such that ipi(x) — > as x — > +oo 
(guaranteed by Statement A.l). Let R(x) G C 2 (x , +c») 7 R(x) — > when x — > +oo and 
R(x),R"(x) G Li(xo,+oo). T/ien i/ie equation 

rf 1 - (Q(x) + R{x))r) = (B.5) 

has unique (up to constant factor) solution r)i(x) such that: 

(i) r)i(x), r][(x) are monotonous functions on some interval [a; +oo) ; a > x ; 

(ii) i]i(x), i][(x) tend to zero when x — > +oo; 

(Hi) the solution r]i(x) when x — > +oo has the asymptotic 

r h (x)=i,i{x)(C l + o{l)) ] r ] ' 1 (x)=rt(x)(C 1 + o(l)) (B.6) 
for some C\ G R. 

Proof: It is straightforward to check that Qi(x) = Q(x) + R(x) also satisfies the 
conditions of Statement A. 2. Then the statement of Lemma B.l follows immediately 
from formula (A. 6). ■ 

Lemma B.2 Let Q(x) satisfies the conditions (Q1)-(Q3) of Theorem 1. Let F{t) G 
C 2 [— S, S] for some S > and satisfies the condition (NL2). Then for any solution ip(x) 
of Eq.(B.l) such that tp(x) — > as x — > +oo there exists some C G R such that the 
asymptotic of ip(x) is given by formula 

X 

e -S(xo,x) r 

= qTJI^ ( C + °(!)) ' S ( x °> x ) = J VWldt x^+oo (B.7) 

x 

The case C = corresponds to zero solution. This asymptotic formula can be 
differentiated at least once. 
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Proof: Due to the condition (NL2) one can introduce f(t) = F(t)/t such that 
/(0) = 0. Suppose that ip(x) is a solution of Eq.(B.l) such that ip(x) — > as x — > +oo. 
Let us rewrite Eq.(B.l) in the form 

r - (Q(x) + R{x))^ = 0, R(x) = -f(4>(x)); 

One can show that the conditions of Lemma B.2 imply the conditions of Lemma B.l 
for R(x). Then the statement of Lemma B.2 follows from Lemma B.l. ■ 

Lemma B.3 Let Q(x) satisfies the conditions (Q1)-(Q3) for some x > Xq. Let 
F(t) G C[— 5, 5] for some 5 > and F(t) satisfies the condition (B.4) for \t\ < 5. Let 
ipi(x) be some fixed nonzero solution of linear problem (A.l), such that ipi(x) — > when 
x — > +oo. Then for any C G R there exists a solution ip+(x; C) of (B.l) such that 

lim ^1 = C (B.8) 

Comment . This statement is converse to Lemma B.2. The assertion of Lemma B.3 
says that for any C there exists a solution of (B.l) which at x — > +oo has the asymptotic 

X 

^ + {x) = (C + o(l)) , S(x , x) = J y/Q({jdt (B.9) 

xo 

Proof: Let ip(x) satisfy the equation (B.l). Introduce the function 

where C is free real parameter and ip\(x) is fixed solution of linear equation (A.l), such 
that i>i(x) — > as x — > +oo. Then \l/(x) satisfies the system of equations 

= $ (B.10) 
$' = -P(x)$-(* + C)./((* + C)^ 1 (a;)), P(x) = ^M (B.H) 

where /(t) = F(t)/t. The existence of solution ^(x) such that ^(x) — > as x — > +oo 
follows from Theorem 9.1, Chapter XII of [19] applied to the system (B.lO)-(B.ll). The 
conditions of this Theorem can be verified straightforwardly ■. 

The following Lemma says that for given C the solution of (B.l) with constant C 
in asymptotic (B.7) is unique. 

Lemma B.4 Let Q(x) satisfies the conditions (Q1)-(Q3) of Theorem 1 for x > xq 
and F(t) satisfies the conditions (NL1)-(NL2). Then for any two solutions, ip a (x) and 
ipb(x), of Eq.(B.l) such that ip a (x) , ipb(x) — > when x — > oo and ip a (x) ^ ipb(x) the 
following inequality holds: 

lim * 1 (B.12) 
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Proof: Assume the contrary, that lim^+oo |^j|y = 1. Then, according to Lemma 
B.2, both solutions ip a (x) and ipb( x ) share the same asymptotics 

e -S(xo,x) e -S(xo,x) 

Ux) = Qv^y {c + o(1)) ' Mx) = gv^y {c + o(1)) ' (R13) 

X 

S(x ,x) = [ y/Qlfydt (B.14) 



Subtracting the second equation from the first one in the system 

< - Q(x)^ b + FM = 0, < - Q(x)*p a + F{i) a ) = 

we conclude that the difference A = ip h — ip a obeys the equation 

A" - (Q(x) + G(ip a (x),ip b (x))) A = 0, (B.15) 

g «^).^))-^»:y (B,e) 

Let us denote R(x) = G(t/j a (x) : ^b(x)). It is tedious but straightforward to show that 
R(x) satisfies the conditions of Lemma B.l (the proof involves Lagrange theorem). Then 
it follows from Lemma B.l that 

X 

A{X) = {Cl + S(x l7 x) = j ^W)dt (B.17) 



XI 



Since A ^0 then C\ ^ and lim^+oo = ^ that proves Lemma B.4. ■ 

Combining the results of Lemmas B.2, B.3, B.4 we arrive at the statements (a) 
and (b) of Theorem 1. The statement (c) follows from the uniqueness of the solution 
t/j + (x;C) for given C (the statement (b)) and the property (NLl)(see e.g. [19]). This 
completes the proof of Theorem 1. ■ 

Appendix C. Proof of Lemmas 1 and 2 from Section 2.4 

Comparison Lemma. Let the functions y(t) and x{t), t G [a; b] be solutions of equations 

yu-g(t,y) = o (C.i) 
x tt -f(t,x) = (C.2) 
correspondingly. Let also the following conditions hold: 

(i) f(t,0> are defined on [a;b] x [A; B] and locally Lipschitz continuous (see 

[24]) with respect to £, £ G (A,B, b may be finite or infinite); 

(n) g(t, > f{t, for any t G [a; 6], £ G [A; B\; 

(Hi) f(t,£) is monotone nondecreasing with respect to £, £ G [A; B}. 

Let A < x(a) < y(a) < B and x t (a) < yt(a). Then x t (t) < y t (t) and x(t) < y(t) 
while A < x(t),y(t) < B or for the whole interval t G [a; b] 
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Proof: Introduce the functions 

f f(t,A), £<A ( g(t,A), £<A 

M0= /(U), A<i<B g(t,t) = l g(t,£), A<i<B 
{ f(t,B), i>B [ g(t,B), £>B 

and let x(t),y(t) be the solutions of Cauchy problem for the equations 

m - g(t,y) = 

x tt - f(t, x) = 

such that x(a) = x(a), Xt(a) = Xt(a), y(a) = y(a), y*(a) = yt(o)- The functions x(t) 
and y(t) coincide with x(t), y(t) when A < x(t),y(t) < B. Consider vectors v, u and 
vector-functions G(i;v), F(t, u) defined as 

For Ui, u 2 G R 2 we define the order relation: Ui < u 2 if each component of u 2 is greater 
or equal than the corresponding component of Ui. The function F(t, u) is monotone 
nondecreasing in the sense that if ui < u 2 then F(t, Ui) < F(t,u 2 ). Then 

= u t - F(t, u) = v t - G(f , v) < v t - F(t, v) 

In this situation Theorem 2 of [23] (see also [24]) states that u(a) < v(a) implies that 
< v(t) for t E [a; b]. This gives the assertion of Comparison Lemma. ■ 

Proof of Lemma 1. Let us denote u(x, C) = d ^' c) . Then u(x,C) satisfies the 
linear equation 

u" + (u-V(x) -3^ 2 (x,C))m = (C.3) 

The conditions of Lemma 1 imply that there is an interval [C+, C\) such that for 
any C G [C\,C\) and any x for which the solution ip + (x,C) is defined, the following 
conditions hold 

ifj + {x;C)>0 (C.4) 
oj-V(x) -3i/jI(x;C) < (C.5) 

Let at least one of the conditions (C.4)-(C.5) breaks at C — C\ i.e. one at least one 
of the signs "<", ">" in (C.4)-(C.5) becomes "=" at some value x — a. It follows 
from asymptotics (15) that there exists b large enough that ip + (x;C) for x > b and 
C G C\\ is monotonic with respect to C in the sense that if C\ > C 2 , Ci ;2 G [C° , C^_] 
and x > b then Ci) ^ i J +( x 'i C2) and V'+( ;;c ; Ci) < V , +( a; i ^2)- This means that for 

a; > 6 and C G [C° , C|] one has u(x; C) > and C) < 0. After changing x — > — rr 
one can apply Comparison Lemma to equations (C.3) and u" = 0. Then u(x; C) > 
for C G [C° , C_^] not only when x > b but everywhere in the interval where u(x; C) is 
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defined. This implies that ' > for C G [C°,C^_] i.e. tp + (a;C) is nondecreasing 

functon of C. Therefore ?/> + (a; C^) < il> + (a;CY) and this is a contradiction to the 
assumption that one of the conditions (C.4)-(C.5) breaks at C — C]_. ■ 

Proof of Lemma 2. Let us apply Comparison Lemma taking 

f(t,o = -^ + e; g(t,t) = -(u>-v(t))t + ? (C.6) 

A = VU; B = +oo (C.7) 

It is straightforward to check that any solution of the equation ip" + Qip — ip 3 = with 
initial data ip(t ) > y/U, ip'(t ) - arbitrary, is singular. Therefore any solution of Eq(20) 
with the same initial conditions is also singular. Due to oddness of the function g(t,£) 
with respect to £ one concludes that any solution of Eq(20) such that ^(*o) < — v^, 
if)' (to)- arbitrary, also is singular. This implies the assertion of Lemma 2. ■ 
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